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ABSTRACT 

Context. Several observational pieces of evidence indicate that specific evolutionary channels which involve Roche lobe overflow 
are not correctly accounted for by the classical Roche model. 

Aims. We generalize the concept of Roche lobe in the presence of extra forces (caused by radiation pressure or pulsations). By 
computing the distortion of the equipotential surfaces, we are able to evaluate the impact of these perturbing forces on the 
stability of Roche-lobe overflow (RLOF). 

Methods. Radiative forces are parameterized through the constant reduction factor that they impose on the gravitational force 
from the radiating star (neglecting any shielding in case of large optical thickness). Forces imparted by pulsations are derived 
from the velocity profile of the wind that they trigger. 

Results. We provide analytical expressions to compute the generalized Roche radius. Depending on the extra force, the 
Roche-lobe radius may either stay unchanged, become smaller, or even become meaningless (in the presence of a radiatively- or 
pulsation-driven wind). There is little impact on the RLOF stability. 
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1. Introduction 

The Roche model has been widely used to infer the out- 
come of binary star evolution. In this model, only the 
gravitational and centrifugal forces are accounted for to 
compute the equipotential surfaces. However if other 
forces are present in the system, like those responsible for 
mass loss, they should be included in the description as 
well. The modification of the usual Roche model was first 
pointed out by Schuerman (1972) in the context of binary 
systems involving early-type main sequence stars with a 
strong wind. The idea was further explored in various 
directions by Kondo & McCluskey (1976), Vanbeveren 
(1977, 1978), Friend & Castor (1982), Djurasevic (1986), 
Zhou & Leung (1988), Huang & Taam (1990), Drechsel 
et al. (1995), Frankowski & Tylenda (2001), Phillips & 
Podsiadlowski (2002) and Owocki (2007). In the present 
paper, we investigate in a more systematic way the sit- 
uations where the Roche model should be modified by 
considering the different physical processes driving the 
stellar winds, that we briefly describe in Sect. 2. 

The Roche model is generally used to answer two 
different questions, namely: (i) What is the flow geom- 
etry? (ii) Is the star filling its Roche lobe? The first ques- 



tion is related to the geometry of the equipotentials (and 
to the Coriolis force). One important modification of the 
equipotential geometry which can arise in the presence 
of an extra force pervading all space (like radiation pres- 
sure), is that the equipotentials open up in the direc- 
tion of the external Lagrangian point, thus possibly al- 
lowing the matter ejected by the mass-losing compo- 
nent to form a circumbinary disc. That issue is impor- 
tant in the framework of the binary evolution involving 
low- and intermediate-mass components where the for- 
mation of a circumbinary disc seems to be very common 
(see de Ruyter et al., 2006; Frankowski & Jorissen, 2007; 
Frankowski, 2009). 

The second question corresponds to the Roche-lobe 
overflow (RLOF) criterion, which involves the compar- 
ison of the Roche radius with the stellar radius. In the 
present paper we show how, depending on the physical 
process driving the wind, the Roche radius may either be 
unchanged with respect to the classical Roche model, be- 
come smaller, or even become meaningless. This will de- 
pend on the value of the extra force at the stellar surface 
(in contrast to the first question, an answer to which re- 
quires the knowledge of the force everywhere within the 
system). In a mass-losing star, the photospheric radius it- 
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self may become ill-defined, thus complicating the use of 
the RLOF criterion (see Sect. 7.4). 

Both issues (the equipotential geometry and the 
RLOF criterion) will be addressed in the present paper 
for the generalised Roche model when radiation pressure 
or pulsations play a role (Sect. 3). The way to correctly 
account for radiation pressure is discussed in Sect. 4. 
Typical values for the radiation pressure at the surface 
of various classes of stars are given in Sect. 5, and the 
corresponding shapes of the equipotentials are displayed 
in Sect. 6. A numerical fit to the Roche radius will be 
provided in Sect. 7.1, generalising Eggleton's formula 
(Eggleton, 1983) to situations where a radiation-pressure 
force is present. The necessity to abandon the Roche-lobe 
concept in the case of stars suffering from radiatively- 
or pulsation-driven wind mass loss is demonstrated in 
Sect. 7.4. Conclusions are drawn in Sect. 8. 

2. The different modes of wind mass loss 

Holzer & MacGregor (1985), Schatzman et al. (1993), 
Lamers (1997), Willson (2000) and Owocki (2004) have 
reviewed the different types of winds existing across 
the Hertzsprung-Russell diagram according to their re- 
spective driving mechanisms. These mechanisms may be 
grouped in three broad classes: radiation-driven winds 
(associated with high-luminosity objects), pulsation- 
initiated winds, and Alfven wave-induced winds. 

The first class includes line-driven winds, operating 
in luminous hot stars (OB stars, and Wolf-Rayet -WR- 
stars as well) on resonance and subordinate lines (Castor 
et al., 1975; Abbott, 1982). Fast winds with terminal ve- 
locities of the order of 500 - 3000 km s~ are generated 
(Kudritzki & Puis, 2000). At the same luminosities as 
O stars, WR stars have larger mass-loss rates, so an- 
other mechanism - like pulsation-driven mass loss, or 
multiple scattering of photons in an optically-thick wind 
with an ionisation stratification - is probably adding to 
radiation pressure on atomic lines (Glatzel et al., 1993; 
Owocki & Gayley, 1999; Nugis & Lamers, 2000; Owocki, 
2004). Dust-driven winds belong to the same category 
as radiation-driven winds, but operate instead in lumi- 
nous cool stars on the asymptotic giant branch and require 
wind densities large enough to couple dust with gas (Gail 
& Sedlmayr, 1987; Lamers, 1997; Wachter et al., 2002; 
Schroder et al., 2003; Sandin & Hofner, 2003; Sandin, 
2008). Wind terminal velocities are small, of the order of 
10-15 km s~ . An important specificity of cool-star winds 
is that the driving force only slightly exceeds the gravita- 
tional attraction, as apparent from their small terminal 
velocities. 

The dust formation requires another process to lift the 
matter high enough above the photosphere, in the region 
where the temperature is less than 1500 K. This can be 
done through shock waves associated with stellar pulsa- 
tion (Bowen, 1988), like in Mira or semi-regular variables. 

Alfven waves are the other main class of waves driv- 
ing winds in stars with open magnetic field lines. This 



latter wind mechanism is important for stars that are not 
luminous enough to have a strong radiation pressure, i.e., 
magnetic A, F, G, K and M stars with luminosities lower than 
about 50 L Q (see Sect. 5). It is the leading candidate to ac- 
count for the solar wind. This driving process has the im- 
portant property to be a local phenomenon which does 
not derive from a potential contrarily to radiation pres- 



sure. 



3. The effective potential 

In the Roche model, the two components of a binary sys- 
tem are considered as point sources in circular orbits and 
in synchronous rotation with the orbital motion. It is then 
possible to define a reference frame in uniform rotation 
about the centre of mass of the system, in which the two 
stars are at rest. When distances are expressed in units of 
the orbital separation, time in units of the orbital period 
and masses in units of the total mass {M\ + Mt), the effec- 
tive potential of the system, including the gravitational 
and the centrifugal potentials, is given by 



O = -£L - 

r\ 
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are respectively the distance of a test particle located at 
(x,y,z) to the primary star (labelled "1" and located at 
x\ — p - 1) and to the companion (labelled "2" and lo- 
cated at X2 - p). The centre of mass is located at the ori- 
gin of the coordinate system. The mass of the test particle 
is supposed to be small enough not to disturb the poten- 
tial. When additional forces deriving from a potential are 
present in the system (other than the gravitational and 
centrifugal ones), the effective potential writes 



O = Op 
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l 



x 2 +y 2 



n 2 

When the extra force is caused by radiation pressure, 
Oextra is easy to evaluate because the radiation force has 
the same l/r 2 dependence on the distance as the gravi- 
tational attraction and it also pervades all space. The va- 
lidity of this simple model is further discussed in Sect. 4. 
The ratio / of the radiation to the gravitational force of 
star 1 

i 



1 dP rad 
f ~ p dr 



(GM, 



(2) 



is independent of r\ provided that the radiation flux at 
position r\ emanating from star 1 - denoted F Vj i(ri) - and 
appearing in the expression 

dP ra d 



dr 



i i r 

ri C Jo 



p F Vj i(n) dv, 
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may be replaced by L v j/(47ir 2 ), so that 

In the above equations, p is the density of the circumstel- 
lar matter, k v is the absorption coefficient per unit mass 
at frequency v, c the speed of light, and L Y \ the lumi- 
nosity of star 1 in the frequency range (v, v+dv). We will 
ignore radiation shielding and assume that the medium 
is optically thin (see also Huang & Taam, 1990; Drechsel 
et al., 1995). If radiation does not reach thermal equi- 
librium with circumstellar matter (i.e., the matter is op- 
tically thin), the radiation luminosity remains approxi- 
mately constant all over the place, and the radiation flux 
can be expressed as F v i(ri) = L V; i/47ir 2 , which may be 
rewritten P Vj i(ri) = P'^PiM) 2 , F* vX being the radiation 
flux emitted per unit surface by the star of radius R\ . If 
that star is tidally distorted, limb- and gravity-darkened, 
a rigourous treatment would imply that Eq. 3 depends 
on the spherical coordinates 9, (p since the flux irradiated 
by star 1 in a given direction defined by 9, (p will now de- 
pend on these variables. The inclusion of these complica- 
tions is beyond the scope of this paper, but they are not 
expected to alter our general conclusions. It is assumed as 
well that k v does not depend strongly on the optical depth 
t in the wind, i.e. that the ionisation state and/or level 
population of the gas do not change. This issue has been 
thoroughly studied in the context of radiatively-driven 
hot-star winds by Abbott (1982), Pauldrach et al. (1986), 
Shimada et al. (1994), Gayley (1995), Lamers (1997) and 
Puis et al. (2000). We thus consider here an idealised situ- 
ation where / is assumed constant. Under these hypothe- 
ses, the effective gravity of the mass-losing star is then 
reduced by a factor (1 - f), the same everywhere, and the 
effective potential becomes 

cD = _ iz£ _ l^f- . (5 ) 

n n 2 

However, even though the equipotentials give an idea 
of the flow geometry, its precise description is only given 
by the equation of motion: 

d 2 r 1 dr 

— = — VP -2fix V<t> + F', (6) 

df 2 p df 

which includes the important Coriolis term. In the above 
expression, P is the gas pressure (this term can be impor- 
tant in case of very dense winds like in Wolf-Rayet stars), 
SI is the angular velocity, F' is the extra force and O is the 
Roche potential (Eq. 1). This equation of motion will not 
be used in the present paper. 

4. How to correctly account for radiation 
pressure? 

The impact of radiation pressure on the shape of the 
equipotential surfaces has been discussed by several au- 
thors (Schuerman, 1972; Huang & Taam, 1990; Drechsel 



et al., 1995; Howarth, 1997; Maeder & Meynet, 2000; 
Howarth & Smith, 2001; Phillips & Podsiadlowski, 2002; 
Owocki, 2007), in slightly different contexts, and there 
has been some controversy as to whether or not this ef- 
fect modifies the Roche geometry. To make the discussion 
clear, one should distinguish three different situations: (i) 
the effect of radiation pressure on the matter outside the 
stars and flowing in the binary system; (ii) the impact of 
the external radiation field (from the companion star) on 
the equilibrium shape of the irradiated star and (iii) the 
effect of the star's own radiation pressure on its equilib- 
rium configuration. 

The first effect is relevant for deriving the flow ge- 
ometry, and the radiation pressure must always be in- 
cluded in the computation. Complications here are only 
related to possible shielding effects. In the present pa- 
per, we consider the situation where radiation pressure 
comes only from the primary and neglect any shield- 
ing by the companion star or by the circumstellar ma- 
terial which remains optically thin everywhere. Case (ii) 
has been considered by Drechsel et al. (1995), Phillips & 
Podsiadlowski (2002) and Owocki (2007) in full details. In 
the present approach, it is not accounted for and would 
require the incorporation of a new (1 — /') factor reducing 
the gravity of the star labelled 2 in Eq. 5. Drechsel et al. 
(1995) and Phillips & Podsiadlowski (2002) have per- 
formed detailed numerical calculations for this situation, 
including geometrical shielding effects (not considered in 
our simple analytical approach) but the impact is usually 
small. Case (iii) is conceptually more intricate and over 
the years has become controversial (Friend & Castor, 
1982; Howarth, 1997; Maeder & Meynet, 2000; Howarth 
& Smith, 2001). To understand why, one should first 
distinguish the optically-thick from the optically-thin 
regimes, following the insightful discussion of Friend & 
Castor (1982). In optically-thick regions, the diffusion ap- 
proximation for the radiative transfer may be used, re- 
sulting in an isotropic pressure with the (thermal equilib- 
rium) value (l/3)«r 4 . Therefore, in the hydrostatic layers 
below the stellar photosphere where this regime holds 
(this excludes WR stars), the effect of radiation should be 
included in the pressure term in the hydrostatic equilib- 
rium equation and the gravitational potential is unmod- 
ified. On the other hand, in an optically-thin region, the 
radiation field is partially decoupled from the gas and 
the stellar flux can be treated in a free-streaming approx- 
imation. In this regime, the radiative flux is radial and 
proportional to 1/r 2 and acts as a repulsive force that re- 
duces the gravitational force by a factor (1 - f). 

The diffusion and streaming approximations hold re- 
spectively deep within the star and in the outer layers 
of wind. Near the photosphere or in the surface lay- 
ers where the hydrostatic equilibrium condition breaks 
down (like in WR stars, see Sect. 7.4) the situation is more 
ambiguous as the von Zeipel's theorem does not apply 
(Howarth, 1997). 

To illustrate this dichotomy, consider Fig. 1 which 
shows the Roche-lobe filling criterion in a giant star at 
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Fig. 1. A giant star with an extended atmosphere (grey re- 
gion, on an exaggerated scale). As discussed in the text, in the 
optically-thin (r < 2/3) layers of the atmosphere, the ratio / 
of the radiative force to gravitational attraction is positive and 
affects the Roche potential (dashed lines). Within the photo- 
sphere, the matter is optically thick, so that / = and the clas- 
sical Roche potential applies (solid lines). The part of the atmo- 
sphere located above its critical Roche lobe will flow through 
the companion (depicted by the arrows). 



(/ = 0) and above (f > 0) the photosphere where the ra- 
diation pressure modifies the potential. The unmodified 
Roche potential (Eq. 1) should be used for the subpho- 
tospheric matter, whereas the modified potential (Eq. 5) 
should be used above the photosphere (grey area on 
Fig. 1). Since the modified potential leads to a smaller 
critical surface (as will be shown in Sect. 7.1), the sit- 
uation may arise where the photosphere does not fill 
its Roche lobe, but the supra-photospheric matter does, 
especially for (super) giant stars with extended atmo- 
spheres. These supra-photospheric shells are removed 
from the star, and in the case of a giant star with a convec- 
tive envelope, this removal will cause the photosphere to 
expand on a thermal time scale (Ritter, 1996). Therefore, 
even though the photosphere itself does not fill its own 
Roche lobe, the consideration of the modified Roche po- 
tential alters the evolution, at least for giant stars. This 
distinction between dwarfs and giants is important, be- 
cause empirical arguments call for the use of a modi- 
fied Roche equipotential for giants, but not for dwarfs. 
Howarth (1997) has given convincing arguments that the 
unmodified potential should be used whenever there is 
evidence for gravity darkening (a direct consequence of 
von Zeipel's theorem), as is the case for dwarf stars in 
binary systems (Anderson & Shu, 1977; Rafert & Twigg, 
1980). The situation is however different for giants: an 
important empirical motivation for considering the case 
/ > is that it may account for the fact that some giant 



stars exhibit ellipsoidal variations (due to non-sphericity) 
despite small (classical) Roche-filling factors. Those el- 
lipsoidal variables would appear quite enigmatic with- 
out the presence of an additional force reducing the ac- 
tual Roche radius and disturbing the stellar shape. As 
it will be shown in Sect. 7.1, the modified Roche radius 
(for / > 0) is smaller than the classical one (for / = 0). 
Therefore, stars nearly filling their Roche lobe with / > 
would be far from filling their lobe if that filling factor 
were estimated with the traditional formula correspond- 
ing to / = 0. This is probably what happens for the 13 
s-type symbiotic systems with ellipsoidal variations de- 
tected by Mikolajewska (2007), which have a (classical) 
Roche filling factor of only 0.4 - 0.5. Based on the anal- 
ysis of the orbital circularisation in a sample of binary 
systems with M-giant primaries, Frankowski et al. (2009) 
find that these giants also do not fill more than ~ 0.5 of 
their classical Roche lobes. Again, the explanation may 
lie in a decrease in the Roche radius below its classical 
value. 



5. Typical / values 

To clarify the ideas, we now evaluate the / factor us- 
ing the Castor- Abbott-Klein (Castor et al., 1975, hereafter 
CAK) theory of line-driven winds. Let g gra v(r) = be 
the gravitational acceleration at distance r from the stel- 
lar centre, and 

fed = Sihd + r(r)) = -^rd + r(T)), (7) 
4n c r L 

be the sum of the radiative accelerations due to Thomson 
scattering (gxh) and to an ensemble of lines (g-rh T(f)), 
where f is the so-called "force multiplier". In the above 
relation, K e = cr e n e /p is the opacity coefficient per unit 
mass for Thomson scattering, cr e is the corresponding 
Thomson cross section (6.65 10~ 25 cm -2 ), n e is the electron 
number density and p is the density of the wind. If the 
wind is fully ionised, K e = <r e ^ (j\ ~ 0.2(1 + X) cm 2 g -1 , 
where mn is the hydrogen mass. The force multiplier 
has been computed by CAK, Abbott (1982), Pauldrach 
et al. (1986), Shimada et al. (1994), Gayley (1995), Lamers 
(1997) and Puis et al. (2000) for various wind thermody- 
namical conditions. It turns out to be very small at optical 
depths close to one because of the "saturation effect": ra- 
diation with frequencies matching those of atomic tran- 
sitions are efficiently removed from the flux at the base 
of the photosphere and no flux at those frequencies is 
left further up in a static atmosphere to accelerate mat- 
ter. Of course, the situation changes in the presence of a 
wind, since the associated velocity gradient will then in- 
duce a Doppler shift of the line frequencies. At low opti- 
cal depths (i.e., higher up in the expanding atmosphere), 
when all lines become optically thin, T(j) becomes very 
large, and reaches asymptotic values up to about 2000 
(see e.g., Table 2 of Abbott, 1982, also Gayley 1995). 
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Table 1. The Eddington parameter T (Eq. 8) for stars along the 
main sequence, adopting K e = 0.35 cm 2 g _1 . Since T < 2000 
(Abbott, 1982), / max = 2000r. Luminosities and masses are from 
Cox (2000). 



Sp. Typ. log (L/L ) M F / ma 



(M Q ) 


05V 


5.95 


60.0 


4.0 10 


-i 


807. 


B0V 


4.77 


17.5 


9.2 10 


-2 


184. 


B5V 


2.97 


5.9 


4.3 10 


-3 


8.68 


B8V 


2.33 


3.8 


1.5 10 


-3 


3.06 


A0V 


1.77 


2.9 


5.5 10 


-4 


1.10 


A5V 


1.19 


2.0 


2.1 10 


-4 


0.42 


F0V 


0.86 


1.6 


1.3 10 


-4 


0.25 


F5V 


0.51 


1.4 


6.4 10 


-5 


0.13 



Table 2. The / parameter for M = 1M Q K and M giants and 
subgiants with C/0= 0.5, at t Ross = 1. 



^ff (K) 


logg 


/ 


5000 





0.97 


4000 





0.43 


3000 





0.18 


5000 


3 


9.6 10~ 4 


4000 


3 


4.3 10~ 4 


3000 


3 


1.2 10- 4 



those stars, the radiation pressure is dominated by the 
contribution of the near-infrared continuum where these 
stars emit most of their radiation. The corresponding / 
values are listed in Table 2. 

For AGB stars with still higher luminosities than 
those considered in Table 2, the radiation driving force 
is now dominated by absorption and scattering in molec- 
ular lines, yielding values of / as large as 0.15 for C-type 
stars with C/O = 2, logg = -1 and T efi = 2500 K (Elitzur 
et al., 1989; Jorgensen & Johnson, 1992). 

In the next sections, the values of / quoted above will 
be applied to different astrophysical situations. It has to 
be made very clear that one should distinguish situa- 
tions involving RLOF (Sect. 7) from situations involving 
modifications of the geometry of the equipotential sur- 
faces far above the photosphere of the mass-losing star 
(Sect. 6). Although in the latter case, the / values pro- 
vided in Table 1 may be used without restrictions, the sit- 
uation is a bit more complicated in the former case. This 
is because the RLOF criterion involves photospheric lay- 
ers, and we have argued in Sect. 4 that non-zero / values 
associated with radiation pressure do not usually alter the 
stellar equilibrium configuration (according to the von 
Zeipel theorem), except in special circumstances involv- 
ing giant stars. The two situations are therefore discussed 
in separate sections below. 



Using the Eddington parameter 



we find that 

/ = r(i + n (9) 

In Table 1, we have listed the / and F values for stars 
along the main sequence. The numbers are obtained by 
considering an optically-thin medium, where T reaches 
its maximum value of about 2000 (Abbott, 1982; Gayley, 
1995). 

The threshold / = 1 is found around spectral type 
AO (corresponding to stars with luminosities less than 
about 50 L Q ), and this result is consistent with empiri- 
cal arguments stating that the peculiarities of Am stars 
cannot survive if the mass loss is too strong (Michaud 
et al., 1983; Michaud & Charland, 1986; Lemke, 1990; 
Babel, 1992). Furthermore, the very existence of Am stars 
demonstrates the importance of the radiative accelera- 
tion on the distribution of various chemical elements 
of these stars (e.g. Hui-Bon-Hoa et al., 2001; Alecian & 
LeBlanc, 2002). 

For late-type stars, the CAK model cannot be used 
because this model does not include the relevant tran- 
sitions, especially molecular lines which may contribute 
significantly to the radiative driving force (Jorgensen & 
Johnson, 1992). Therefore, radiative accelerations for late- 
type stars have been taken from the MARCS model atmo- 
spheres (Gustafsson et al., 2008) for K and M giants. For 



6. The modified Roche equipotentials 

Fig. 2 presents cross-cuts through the effective equipoten- 
tial surfaces (Eq. 5) along the line joining the two stars, for 
H = 1/3, i.e. the mass-losing star (located at x\ = -2/3) is 
the less massive component, and for three different val- 
ues of /: / = (solid line), / = 0.7 (short-dashed line), and 
/ = 1.1 (long-dashed line). The first case corresponds to 
a situation with no extra-force, the second to a situation 
where the extra-force is present but not sufficient to drive 
mass loss (note that the Roche lobe is becoming smaller 
as compared to the case / = 0), and the third case cor- 
responds to a situation where the extra-force is actually 
driving the mass loss. Note that in this case, there is no 
Roche lobe any longer. 

Fig. 3 presents the different families of equipotentials 
as a function of / and when / < 1. Panels b and c corre- 
spond to the critical configurations where the Lagrangian 
points L\ and L3 or L2 and L3 are located on the same 
equipotentials. These situations are encountered for the 
specific values of / denoted fi(fi) and fii/j), respectively. 
These functions are approximated by the following ex- 
pressions with a relative error smaller than 1 % for f\ (jj.) 
and 2% (when fi < 0.97) for f 2 Qi): 

fiUi) = (1 -2yu)(l -0.80yL/ I/2 + 0.82ju) (10) 
f 2 (M) = (1 - H) (l - 0.72yu 1/2 + 0.66|/ 3/2 + n u ) . (11) 

The term /j. 11 is necessary to obtain a good fit to fiin) 
when fi approaches 1 . 
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Fig. 3. Section in the orbital plane of the Roche equipotentials for a binary system with a dimensionless mass /j = 1/3 and / = 
(panel a), the critical values f\ = 0.27 (panel b) and f 2 = 0.48 (panel c), obtained from numerical integration. The case / = 0.6 is 
presented on panel d. The primary star, located at x = fi - 1, is depicted by the cross. The Lagrangian points are labelled 1, 2 and 3. 
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X 

Fig. 2. Cross-sections through the effective equipotential sur- 
faces (Eq. 5) along the line joining the two stars, for /j = 1/3. 
The mass-losing star (located at x\ = -2/3) is the less massive 
component and has a (dimensionless) radius of 0.1. Three cases 
are depicted: / = (solid line); / = 0.7 (short-dashed line; note 
that the Roche lobe is becoming smaller as compared to the case 
/ = 0); / = 1.1 (long-dashed line: there is no Roche lobe any 
longer). 



The different regions delineated by f\ and fa in Fig. 4 
correspond to different equipotential topologies. Systems 
with / < fu f = fu /i < f ^ h and / > fz are topo- 
logically similar to cases displayed in Fig. 3 in panels 
a, b, c and d, respectively. In the last case (panel d), the 
Roche lobes of the two components do not belong to the 
same equipotential. In contrast to the standard case, the 
matter ejected by the primary is not necessarily trans- 
ferred directly into the Roche lobe of the companion, but 
all or a fraction of it may instead feed a circumbinary 
disc. As mentioned in Sect. 1, this issue is important as 
such discs are very common in binaries involving low- 
and intermediate-mass components. The Coriolis force, 
which is not conservative, plays also a major role in the 
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Fig. 4. The functions flip) and f 2 {n) (see text), obtained by nu- 
merical integration (dots) and their analytical fit (solid curves) 
as expressed by Eqs. 10 and 11. The three regions delineated by 
the curves in the (f,fi) plane correspond to three different topo- 
logical structures for the equipotentials (see text and Fig. 3). 
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Fig. 5. The normalised Roche radius Rr/ci (where a is the orbital 
separation) as a function of the mass ratio (q = Mj/M 2 ) for the 
cases / = (solid line), / = 0.5 (long-dashed line) and / = 0.9 
(short-dashed line). 



formation of such discs, as shown by numerical simula- 
tions (Theuns & Jorissen, 1993; Mastrodemos & Morris, 
1998, 1999; Sytov et al., 2009). 
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7. The RLOF criterion 

7. 1. The modified Roche radius 

A generalisation to the cases < / < 1 of the Eggleton 
(1983) formula for the radius of a sphere with the same 
volume as the Roche lobe, is given by 



K(q,f) = 



R R (qJ) 



Mf) q 



2/3 



(12) 



a B( / ') q^ + In (1 +C(/V / 3 ) 
where 

A(f) = (1 - /) 1/3 (0.49 + 0.25/ + 0.35/ 2 

-0.59/ 3 + 0.37/ 4 ), 
B(f) = 0.6 + 0.3/, 
C(/)= 1+/, 

q = M1/M2 is the mass ratio of the two components 
(star 1 being the one with non-negligible radiation pres- 
sure) and a is the orbital separation. For / = 0, the 
Eggleton (1983) expression is recovered. Fig. 5 shows the 
Roche radius as a function of the mass ratio q for the val- 
ues / = 0, / = 0.5 and / = 0.9. The Roche radius is clearly 
smaller in the presence of radiation pressure and is re- 
duced by a factor of ~ 2 between the cases / = and 
/ = 0.9. The Roche lobe vanishes (and its radius thus 
goes to zero) when the radiation force becomes equal to 
the gravitational attraction (i.e. / = 1). 

The analytical fit to the Roche radius as approximated 
by Eq. 12 is accurate to better than 3% over the extended 
range 0.1 < q < 00 and < / < 0.95. Outside this parame- 
ter range, the relative error is less than 7%. The fit is based 
on numerical results obtained with the method outlined 
by Huang & Taam (1990). We emphasize once again that, 
according to the discussion of Sect. 4, the Roche radius 
expressed by Eq. 12 applies only to stars with an ex- 
tended atmosphere where the free-streaming approxima- 
tion holds for the radiation field. 



On the other hand, thermal equilibrium is recovered on 
the Kelvin-Helmoltz time scale (t kh ) and is characterised 
by the exponent The adiabatic and thermal stability 
conditions become respectively 



ft < ftd 

and 

ft < ft> • 



(15) 



(16) 



If neither of these conditions is satisfied, mass transfer 
proceeds on the fastest of the two timescales. 

To estimate £r, Soberman et al. (1997) (see also 
Jorissen, 2003) assume that a fraction a of the mass lost in 
the wind escapes to infinity, that a fraction q is accreted 
by the companion and the remaining 5 — 1 - a — 27 goes 
into feeding a circumbinary disc of radius a r - y 2 a. In this 
framework the Roche-lobe mass-radius exponent writes 



ft = 2 



1 + q 



+ y5(\ + q) -A + 



q 



1 + q 



+2(qq 



ding 



(17) 



where A ~ 0.05 according to the hydrodynamical simula- 
tions of Theuns et al. (1996). 

When an extra force is present, the term (dln'R/dlng) 
increases and may possibly destabilise the system. 
However, the effect is small (< 10% for / ranging be- 
tween and 0.9) and the stability condition remains al- 
ways dominated by the first three terms of Eq. 17, which 
depend on the efficiency of mass transfer through the 
parameters q, 5 and y. This is apparent from Fig. 4 of 
Soberman et al. (1997) or Fig. 9.14 of Jorissen (2003). The 
direct effect of the extra force on RLOF stability is conse- 
quently negligible. 



7.2. RLOF stability 

We now evaluate the impact of the above modifications 
on the RLOF stability. In the Roche model, the stability 
condition imposes that, when the star fills its Roche lobe, 
subsequent mass loss does not lead to a runaway situa- 
tion. It is expressed by the condition £ R < f , where 



ft s 



A\nR R 
din Mi 



and 



dln/?i 
din Mi 



(13) 



(14) 



are the Roche-lobe mass-radius exponent and the mass- 
radius exponent of the donor, respectively. A star re- 
sponds to mass loss on two timescales. The immediate 
response is on the adiabatic time scale (T dyn ), after which 
hydrostatic equilibrium is restored but negligible heat 
transport has occurred. The mass-radius exponent char- 
acterizing this adiabatic readjustment is denoted by £ a( j. 



7.3. Critical period 

An interesting application of Eq. 12 is related to the crit- 
ical period f cr ; t , the orbital period below which RLOF oc- 
curs, given by (Eggleton, 2006) 



At?R- 
GM 



0.1159 



3\l/2 



M(M ) 



1/2 



K(q,f)- 3 ' 2 (days), 



(18) 



where M = Mi + M%. f CI i t increases strongly when an ex- 
tra force is present as shown in Fig. 6. This means that 
RLOF, either stable or not, would occur at a longer or- 
bital period than inferred from the size of the classical 
Roche lobe. Such an effect of the extra force (with an esti- 
mated / = 0.65 - 0.85) shrinking the effective Roche lobe 
was offered by Frankowski & Tylenda (2001) as an ex- 
planation for the symbiotic stars "avoiding" to fill their 
Roche lobes (Miirset & Schmid, 1999). As noted in Sect. 4, 
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Fig. 6. Critical period P clil in days from Eq. 18 as a function of 
the stellar radius R(R Q ) for / = (solid line), / = 0.5 (long- 
dashed line) and / = 0.9 (short-dashed line). In each case, Mi = 
1.2 M s and M 2 = 0.6 M G (i.e., q = 2 or ju = 2/3). 

it was also discovered that a number of symbiotics ex- 
hibit ellipsoidal variations, apparently being tidally dis- 
torted despite not quite filling their classical Roche lobes 
(Mikolajewska, 2007). Observations of the orbital circu- 
larisation in a sample of binaries with M giant primaries 
supply another possible example of this effect: these gi- 
ants do not fill more than ~ 0.5 of their classical Roche 
lobe (Frankowski et al., 2009). If this interpretation of 
these systems properties being modified by the extra 
force is correct, the / > 0.7 values implied for these cool 
giants are significantly higher than the estimates given 
for absorption/scattering in molecular lines in the pho- 
tospheres of AGB stars (see Sect. 5). They are closer to 
the effective / ~ 1 values inferred from the pulsation- 
driven winds of AGB Mira stars. This reflects the fact 
that radiation pressure is not the only energy and mo- 
mentum source available in the photospheres of giant 
stars (e.g., even non-pulsating cool giants exhibit bulk 
motions, manifesting as radial-velocity jitter; Gunn & 
Griffin, 1979). 




-1.5 -1 -0.5 0.5 1 

X 



Fig. 7. Roche potential surface (upper panel) and equipotentials 
on the orbital plane (bottom panel) for /j = 2/3 and / = 1.1. The 
mass-losing star, located at x = fi - 1 and depicted by the cross, 
has a radius of 0.1 (in units of the orbital separation). The only 
Lagrangian point is L 2 . 



7.4. No RLOFfor / > 1 

We showed in Fig. 2 that when / > 1, there is no longer a 
critical Roche surface around the mass-losing star, so that 
the very concept of RLOF becomes meaningless. 

For a star where, e.g., the radiation pressure is large 
enough to expel stellar material, the potential corre- 
sponds to a net repulsive force (long-dashed-line in 
Fig. 2) and the Roche lobe around the mass-losing star 
has no meaning any longer! Similarly the stellar radius 
needs to be re-defined, especially in the case of optically- 
thick winds (like for WR stars), when the photospheric 
radius (corresponding to an optical depth r = 2/3) falls 
within the wind (de Loore et al., 1982; Baschek et al., 1991; 
Moffat & Marchenko, 1996). This property is accounted 
for in recent stellar models which correct the stellar ra- 
dius using extrapolation of the wind expansion law in 
the optically-thick region (Langer, 1989; Hamann, 1993). 



A further consequence of the absence of a Roche lobe 
around radiatively driven mass-losing stars is that there 
is no dramatic change in the mass-loss regime from wind 
mass loss to RLOF, as the latter is now ill-defined. In 
fact, several authors have already promoted this idea of a 
smooth transition of the mass loss rate from the wind to 
the RLOF regime (Tout & Eggleton, 1988; Frankowski & 
Tylenda, 2001). 

An important issue related to mass transfer is to eval- 
uate whether or not a common envelope will form and 
this outcome depends on the geometry of the equipoten- 
tials. Contrarily to the situation prevailing during clas- 
sical RLOF, the mass lost by the wind is not necessarily 
injected into the Roche lobe of the companion as illus- 
trated in Fig. 7. In particular, a substantial fraction of the 
wind can avoid the companion's Roche lobe and instead 
be used to form a circumbinary disc. 
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Fig. 8. Wind velocities v wind (ri) (solid, short- and long-dashed 
lines and left-hand scale) and escape velocity (dot-dashed line) 
from model calculations of Bowen (1988) for a Mira star of M = 
1.2 M e , R = 270 R B and L = 5415 L Q , including pulsation and 
radiation pressure on dust. The function /(n ) (dotted line and 
right-hand scale) is derived from Eq. 20. 



7.5. Pulsation-driven winds: The case of Mira stars 

The process driving the wind in Mira stars seems to 
start with pulsation-induced shock waves that lift the 
matter high enough above the photosphere for dust to 
form. Since all photospheric particles will feel this up- 
ward force, it should be included in the effective poten- 
tial. However, in the case of momentum transfer from 
shock waves, there is no simple mathematical expres- 
sion for the extra acceleration g ext , unlike the case of the 
radiation-pressure force. Nevertheless, it is possible to in- 
fer the run of g ext as a function of the distance r\ from the 
stellar surface using model predictions for the wind ve- 
locity (see e.g., Bowen, 1988; Willson, 2000). In the steady 
state approximation, the acceleration dv/df for a wind 
particle writes 



, grav 



dv 

dv/df = v — 
dri 



(19) 



The above relation implicitly assumes that the wind ve- 
locity function is derivable. It thus requires to smooth the 
discontinuities associated with the shock waves (thick 
line in Fig. 8). By suppressing the non-conservative char- 
acter of the shock waves, this case becomes similar to that 
of a conservative force deriving from a potential <J> e xt- The 
fact that Mira stars have relatively low-velocity winds 
implies that this extra force is mostly used to lift the mat- 
ter out of the potential well. The run of / with distance 
then writes 



fin) = 



gext(n) 
ggravCn) 



= 1 + 



dv 



(20) 



Because Mira stars are pulsating, their wind-velocity 
curves are time-dependent (Bowen, 1988). In order to 
mimic this time variability, high-frequency spatial and 
temporal variations are added to the smooth, long-range 
velocity curve. The high-frequency component corre- 
sponds to a sinusoidal curve whose amplitude decreases 
exponentially with distance from the stellar surface. This 



$ _2 




Fig. 9. Same as Fig. 2 in the case of a Mira-type wind (solid, 
short- and long-dashed lines), derived from the wind velocity 
curves displayed on Fig. 8 and using Eq. 20. The star is assumed 
to have a radius of 0.1 (in units of the orbital separation), and 
fj = 2/3. The dotted and dot-dashed lines correspond respec- 
tively to the cases / = 1 and / = 0. 



high-frequency component is moreover phase-shifted by 
7r/2 (solid line in Fig. 8) and n (long-dashed line) to mim- 
ick the temporal evolution (shock-wave propagation). 
The corresponding potentials are represented in Fig. 9 
(solid, long- and short-dashed lines). Inside the stellar ra- 
dius, matter is supposed to be in hydrostatic equilibrium, 
with gas and radiation pressure balancing gravitation, so 
that matter is at rest on average. At different times of the 
shock-wave propagation, the potential can be repulsive 
(solid and long-dashed curves), allowing for the ejection 
of matter, or becomes attractive (short-dashed line), tem- 
porarily preventing mass ejection. 

In our calculations, / remains on average only slightly 
larger than unity which confirms the fact that the driv- 
ing mechanism is mainly used to work against the grav- 
itational attraction. Finally, for binary systems involv- 
ing Mira stars, the companion does influence the mass- 
losing star by altering the physical processes driving its 
wind. Frankowski & Tylenda (2001) for instance stress 
that the empirical formulae fitting red-giant mass-loss 
rates (Reimers, 1975; Arndt et al., 1997) depend on the 
surface gravity, and consider how gravity of the mass- 
losing star will be altered, both because of the straightfor- 
ward addition of the companion gravitational attraction 
and because of the tidal distortion. Furthermore, in Mira 
stars, the gravitational field of the companion may also 
alter the source of the mass-loss, by disturbing the forma- 
tion and the properties of the shock waves. All these as- 
pects contribute to making the determination of / rather 
uncertain. 
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8. Conclusions 

This paper presents and analyses the Roche potential 
modified by the presence of an extra force associated with 
radiation pressure or pulsation. The magnitude of this 
perturbing force is quantified by the parameter / which 
represents the ratio of the extra-force and the gravita- 
tional attraction (Schuerman, 1972). An estimate of this 
parameter for main sequence, RGB, AGB and Mira stars 
is also provided. 

For < / < 1, the Roche potential may be substan- 
tially modified. In particular, if / > f\ the deformation of 
the equipotentials allow the matter ejected by the mass- 
losing star to go into a circumbinary disc. As the extra 
force (/) becomes stronger, the Roche radius decreases, 
favouring RLOF mass transfer. 

Numerical fits and generalisation of the Roche radius 
are provided in this paper for f < 1 . It is shown that the 
effects of the extra force on the RLOF stability is negligi- 
ble. 

For / > 1 , the Roche lobe has no meaning any longer. 
Such situations occur in luminous stars where radiation 
drives the mass loss or in pulsating giant stars. In this 
latter case, the recurrent deposition of momentum by the 
shock waves in the atmosphere allows matter at the sur- 
face to be expelled. 

The consideration of a modified Roche lobe is (di- 
rectly or indirectly) supported by numerous observa- 
tions (like the frequent occurrence of circumbinary discs 
in post-mass-transfer systems, and the small classical 
Roche-filling factors derived for symbiotic or M giants 
despite their ellipsoidal variability or circular orbit...) and 
should be taken into account. 
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